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Noncontact dissipation reveals critical central peak in SrTiOa 
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The critical fluctuations at second order structural transitions in a bulk crystal may affect the 
dissipation of mechanical probes even if completely external to the crystal surface. Here we show 
that noncontact force microscope dissipation bears clear evidence of the antiferrodistortive phase 
transition of SrTiOs, known for a long time to exhibit a unique, extremely narrow neutron scattering 
“central peak”. The noncontact geometry suggests a central peak linear response coupling connected 
with strain. The detailed temperature dependence reveals for the first time the intrinsic central peak 
width of order 80 kHz, two orders of magnitude below the established neutron upper bound. 


Second order structural phase transitions leave a clear 
mark in all thermodynamical, mechanical, equilibrium 
and non-equilibrium properties of bulk crystals. It was 
proposed some time ago that the critical fluctuations 
should also leave a footprint in the frictional dissipa¬ 
tion of external mechanical probes such as an atomic 
force microscope (AFM) when temperature crosses the 
phase transition in the underlying bulk [T]. The recent 
successful detection of a superconducting transition in 
the linear response mechanical dissipation of a noncon¬ 
tact, pendulum-type AFM tip hovering more than one 
nm above the sample surface [2] suggests that contin¬ 
uous structural transitions might also be detectable in 
this manner. Here we present a first realization of this 
idea, with direct application to a most classic example, 
the antiferrodistortive transition of SrTiOs just above 
100 K. At this phase transition the high temperature 
ideal cubic perovskite crystal structure becomes unstable 
against a zone-boundary phonon-like displacement of the 
ions, leading to a cell doubling and a tetragonal I4/mcm 
symmetry at lower temperatures. This exquisitely sec¬ 
ond order “displacive” transition historically provided a 
clean realization of nonclassical critical exponents |3]. A 
very intriguing feature of this system, originally uncov¬ 
ered by neutron scattering, and later confirmed by other 
techniques, is the so-called “central peak” [IHZ]. Very 
close to the critical transition temperature Tc, inelastic 
neutron spectra showed, besides ordinary critical fluc¬ 
tuations - which proliferate and soften but never reach 
zero frequencies - a strikingly narrow peak (less than 
the 6 MHz width resolution) centered at zero frequency, 
whence the name. The central peak (CP) intensity ap¬ 
peared to obey the static critical exponents of the tran¬ 
sition, but despite considerable efforts the actual nature 
and width of the central peak were not uncontroversially 


established [S]. 

Here we show that noncontact pendulum AFM dissipa¬ 
tion, measured far from actual contact with the surface, 
reveals for the first time a structural phase transition, 
and it does so by revealing the CP of SrTiOs. A linear 
response analysis shows that the CP-related mechanical 
loss peak is as narrow as 80 kHz, a frequency orders of 
magnitude below the neutron established upper bound. 
Moreover, even if it cannot strictly determine the inti¬ 
mate nature of the CP, the mechanical coupling suggests 
a connection with critical fluctuations of strain, which 
are known to be associated with those of the main anti¬ 
ferrodistortive order parameter [2]. 

The probe consisted of a very soft, highly doped silicon 
cantilever (ATEC-Cont from Nanosensors) with spring 
constant k = 0.1 N/m, suspended perpendicularly to the 
surface with an accuracy of 1° and operated in the so- 
called pendulum geometry where the tip vibration de¬ 
scribes an arc parallel to the sample surface. The pecu¬ 
liarity of this technique is to detect phenomena, in this 
case phase transitions, which happen in bulk, by means 
of a non-invasive, ultra sensitive and local surface probe, 
as opposed to traditional probes such as neutrons and X- 
rays which invade the bulk in a much more global fashion. 
Moreover, the pendulum AFM is a kilohertz probe, sen¬ 
sitive to phenomena and to fluctuations that may take 
place on a much slower time scale than that accessible 
with neutrons or X-rays. The oscillation amplitude A 
of the tip was kept constant to approximately 5nm us¬ 
ing a phase-locked loop feedback circuit. The cantilever 
was annealed in UHV up to 700° C for 12 h, which re¬ 
sults in removal of water layers and other contaminants 
from both the cantilever and the tip. After annealing the 
cantilever quality factor, frequency and internal damping 
were equal to Q = 7- 10^, /o = 11 kHz and Wq = 2 • 10“^^ 
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kg/s, respectively. The annealing is also known to reduce 
all localized charges on the probing silicon tip [2] , which is 
neutral, since the tip-sample contact potential difference 
was compensated (V=Vcpd) during the experiment. 

Fig. (a) shows the power W{T) dissipated by the 
pendulum AFM as a function of temperature at different 
spots on the SrTiOa surface and at different tip-sample 
distances, as measured by the shift of the resonance fre¬ 
quency A/. The dissipation is inferred from the stan¬ 
dard expression mw = w^ (Axc(^)/Axc.O - /(^)//o) 
in terms of the measured distance-dependent excitation 
amplitude Aexc(-z) and resonance frequency f[z) (where 
f{z) = fo + A/) of the cantilever, the suffix 0 refer¬ 
ring to the free cantilever. Since the Young modulus 
of the silicon cantilever is temperature dependent also 
the frequency of the free cantilever changes as a func¬ 
tion of temperature - A/o(T) [TT]. In a temperature de¬ 
pendent experiment the total change of the frequency is 
Af{T,z) = A/o(T) -f Af{z) where Af{z) is the (neg¬ 
ative) frequency shift due to the tip-sample interaction. 
The tip-sample distance z was accurately controlled by 
means of feedback loop regulating the z position in such 
a way that A/ was kept constant, (see Supplementary 
Material). 

Data at the large distance z = 12 nm, corresponding to 
A/ = —10 Hz, show a dissipation peak which is barely 
visible, corresponding to an exceedingly weak van der 
Waals tip-surface interaction. All other spectra, taken at 
closer distances, exhibit a narrow dissipation peak at a 
temperature between 114 and 118 K depending upon the 
surface spot investigated, reflecting local changes of Tc 
determined by inhomogeneous heavy Nb doping, surface 
oxygen vacancies, and/or stress irregularities. The 105K 
transition temperature of stress free pristine SrTiOa is 
notoriously shifted by Nb doping and formation of oxygen 
vacancies [12) . At surfaces, moreover, Tc may under suit¬ 
able conditions show differences of tens of degrees with 
respect to the bulk, as seen on SrTiO 3 ( 110 ) [T3]. The 
dissipation peak in this raw data provides a first quali¬ 
tative confirmation of the suggested connection between 
critical structural fluctuations and AFM dissipation [1]. 

Fig a (b) shows a scanning tunneling microscope 
(STM) atomically resolved surface topography of the 1% 
Nb doped SrTiO3(001) surface taken at low tempera¬ 
ture r = 5 K. The flat terraces are obtained after a 2h 
annealing to 1000° C under ultra high vacuum (UHV) 
m US]. Detailed STM images (see also Supplementary 
Material) show dark spots (surface defects, perhaps O va¬ 
cancies m) and bright features, decorating what could 
be edge dislocations m or other domain walls. 

We now consider the origin of the pendulum AFM loss 
process. For a start, the tip is sufficiently far from the 
surface to guarantee that only van der Waals (vdW) (or 
electrostatic if charges were present) tip-substrate inter¬ 
actions are relevant. Pure SrTiOs is an insulator and the 
coupling of a neutral tip must be phononic |17j . Re- 


(a) 




FIG. 1. (a) - Experimental AFM dissipation IT as a func¬ 

tion of temperature. Raw data, taken at different surface 
spots and different tip sample distances z. The sharp peak 
corresponds to the critical temperature of SrTiOa in the bulk 
region under the tip. (b) - Low temperature (T=5K) STM 
image of SrTiO3(100) surface. The image is obtained at con¬ 
stant current I=10pA and bias voltage U=1V. More details in 
Supplementary) (c) - The distance dependence of dissipation 
W, taken as the maximum of the peak shown in Fig.1(a), at 
four different spots on the sample. A fit to the experimental 
data W oc z~^ is shown in red, with p ~ 4.2. This exponent 
is close to the value p = 4 expected for phononic dissipation, 
as appropriate for coupling to acoustical surface fluctuations 
of an insulating bulk material. 


sistivity measurements of 1%-Nb doped crystals exhibit 
conducting behaviour, however with a carrier density of 
about 10^°cm“^ [TH], orders of magnitude below that of 
a good metal. 

Moreover, Auger electron spectroscopy on SrTi 03 - 
(2x2) surface has suggested that the Nb presence is negli¬ 
gible in the near-surface region [Tn| , so that the low level 
metallicity due to Nb doping can be considered irrele¬ 
vant in our experiment. Fig. (c) shows the maximum 
dissipation value against tip-sample separation. For a 
spherical tip oscillating above a solid surface the dissi¬ 
pation is proportional to A^(z), where F{z) is the static 
force resulting from tip-sample interaction. The vdW 
interaction yields a static force F{z) oc z“^, so that the 
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dissipation due to creation of phonons in the solid (acous¬ 
tic phonons in this case, corresponding to the oscillating 
strain wave under the tip sketched in the inset in Fig. 
should vary as z' HE]. The experimental distance 
dependence, is indeed best fit by z~p with p ^ 4.2 is in 
excellent agreement with that expectation. 



the term hujnB{u>,T), with ub the Bose function, in 
the experimentally relevant regime huj kBT. Us¬ 
ing the form by Shapiro et al. |5] which accu¬ 
rately describes neutron scattering, the order param¬ 
eter (zone boundary) susceptibility can be written as 

xiq,u},T) = [U^(q) — -I-n(( 7 ,a;, T)] where U is 

a bare soft phonon frequency far from the transition 
and n ^ A(T) — iujTo{T) is a self energy renormal¬ 
ization from anharmonic effects (we shall from now on 
drop the wave-vector q dependence of these quantities). 
This simple form of 11 would lead, in the standard text¬ 
book description of a displacive transition to a T- 
dependent shift of U, resulting in a lorentzian peak in 
Imx(w) at uJooiT) = + A(T), of width Tq, such 

that ujoo(T) —>■ 0 at T = T^. However, the neutron data 
of SrTiOa showed that phonon softening is incomplete, 
Woo(2c) « 0.5 meV, but accompanied by an extra fea¬ 
ture centered at some very-low-energy Wiow) the central 
peak, phenomenologically captured [5] by an additional 
contribution to the self-energy H 

n(w, T) = A(T) - *a;ro(T) - . (2) 

1 - Zw/Wiow 


FIG. 2. Imxcp(w), the low-energy central peak component 
of Imx(tj) in log-log scale for various temperatures, showing a 
peak at tjpeak(T) ~ that moves towards 0 as T —> Tc. (Up¬ 
per inset) Sketch of the full Imx(a)) in linear ca-scale, showing 
the broad soft-phonon lorentzian at ujoo with the sharp low- 
energy central peak. For clarity we used here a high value 
of caiow to show both peaks on the same scale. (Lower inset) 
Cross section of 3D simulation of a tip perturbing a semi¬ 
infinite crystal through a vdW interaction. The tip (red dots) 
is shown as a truncated pyramid where every atom exerts a 
—C/r® vdW potential on crystal atoms (blue dots), that are 
held together by a harmonic potential. Arrows (magnified for 
clarity) represent on a log scale the atom displacements from 
the relaxed positions. 


For w ~ LUao ^ Wiow one recovers the usual soft-phonon 
lorentzian peak at uJoo, but for uj < wiow a second peak 
appears, well approximated by (see upper inset of Fig.[^ 

^ ^ ^ LUlo^6‘^{T) UJ 

ujtiT) a;2 + Kww2(r)/a;^(r)]2 ’ 

( 3 ) 

where ujq{T) = uj1^{T) — 6‘^{T) is the quantity that actu¬ 
ally vanishes as T — Tc- Indeed, the static susceptibility 
can be shown to be simply related to ojq 



duj Imx(w) 

TT UJ 


1 




( 4 ) 


We can now directly relate the observed dissipation to 
the the critical central peak of SrTiOa. The noncontact 
tip vibrating at / « 11 kHz and at large distance is a very 
weak perturbation on the underlying SrTiOa. Thus, we 
can make use, rarely appropriate in nanofriction, of linear 
response theory. Moreover since the AFM perturbation 
affects a sufficiently large portion of SrTiOa, we can ap¬ 
proximate its response by means of the bulk response of 
the material HEj- The dissipated tip energy per cycle 
is, in linear response, proportional to the imaginary part 
of the bulk lattice susceptibility Xi in Ih® form 

Wiuj,T) = Wo + akBTlmxiuj,T) , ( 1 ) 

where Wq is the dissipation of the free cantilever (T- 
independent in the considered temperature range), x is 
an appropriate momentum average of the lattice suscep¬ 
tibility x(<ZiW,T), a is a positive, distance-dependent 
constant and the temperature factor originates from 


The divergence of the order-parameter susceptibility x 
with an exponent 7 , as the reduced temperature t = 
\T — Tc\/Tc goes to 0 is a standard result of the the¬ 
ory of critical phenomena. The critical behavior of 
SrTiOa is in the 3D-Ising universality class, for which 
7 ~ 1.24 [51]. The low-energy susceptibility Imxcp(w) of 
Eq. (H displays a sharp peak at a frequency a;peak(T) = 
wiow<^(T)/w^(T), which moves towards 0 as T — >• Tc- 
We can now consider the temperature dependence of the 
linear response AFM dissipation at the fixed and very 
low oscillation frequency Wexp = Stt/. As T ^ Tc from 
above the dissipation will increase, roughly as , be¬ 
cause Wexp ^ Wpeak(’E) to reach saturation value 

at T = r such that Wpeak(’E) ~ Wexp- Essentially T 
(here about 1 K above Tc) is the temperature below 
which CP fluctuations average out. Correspondingly, be¬ 
low T the dissipation levels off as we can essentially take 
Imxcp(wexp) « a;iow<5^(T)/[a;expW^(P)], which depends 
very mildly on T. (experimental values for S^(T) and 
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uJoo{T) given by j^). We finally obtain an overall pre¬ 
dicted critical form for the AFM dissipation: 


W = Wo + 


U 

1 -b ’ 


( 5 ) 


where U and V are positive constants. (In the notation 
of Ref. [5], = 7^Wg(r)/(a;^(T)a;gxp) and at low t 

the relevant dependence on temperature is given by the 
a;Q(T) term. ) 



FIG. 3. Experimental dissipation W above Tc as a function 
of temperature (black dots). Inset: the same data in a lin¬ 
ear temperature scale, showing data both below and above 
Tc. Red and green lines: fit above and below Tc according 
to Eq. pink dashed line: fit above Tc excluding the high 
plateau data. The value of oJiow is 156kHz (solid curve) and 
44kHz (dashed curve), both within a factor 2 of 83 KHz, the 
value obtained from a simple ht of the saturation temperature 
(see text), t is roughly the temperature where, upon cooling 
from above Tc the dissipation levels off to a plateau. Hori¬ 
zontal error bars corresponding to St ~ 10“"^ are small and 
omitted. 

Fig. [3] shows on a log-log scale the data for W — Wo 
at A/ = -40Hz {z - 3.5 nm) for T > ~ 117.58 K. 

Considering the experimental uncertainty mainly due to 
noise in the dissipation signal, a slope provides a 
good fit well above Tc, followed by a saturation when 
t <i ^ 10“^. Taking from Ref [5] « 0.3 meV^ and 

ujQ{t) « 0.04 meV^, we finally observe that this satura¬ 
tion of AFM dissipation determines the low-energy width 
parameter wiow as wiow = (t)/wg (t) ~ 83 kHz. 

We draw in summary four conclusions. First, bulk 
structural phase transitions are indeed revealed by AFM 
dissipation, as was predicted. [T] Strikingly, in the present 
noncontact realization, this is realized without literally 
touching the crystal. Second, the pendulum AFM dissi¬ 
pation picks up precisely the long debated central peak 
fluctuations, here responsible for the dissipation at the 
extremely low AFM pendulum frequency of 11 kHz. 


Third, the unknown breadth wiow of the central peak in 
the dynamical structure factor S{uj) = Imx(a;)/a; now 
obtained as an intrinsic property of SrTiOa is about 
80 kHz, well below the upper bound set by the neu¬ 
tron resolution limit of 6 MHz. This CP width is man¬ 
ifested in AFM dissipation as a peak at a;peak(T) = 
wiow<j-'o(T)/a;^(T) « 3.2 MHz. Fourth, the noncon¬ 
tact, large distance tip-surface coupling elicits a phononic 
dissipation attributable in turn to a slowly varying tip- 
induced strain, and not to the primary antiferrodistortive 
order parameter, to which the far away tip and its motion 
cannot directly couple. While this realization does not 
reveal by itself the intimate nature of the CP, which re¬ 
mains open to discussion [8] it does show that the exceed¬ 
ingly slow critical CP fluctuations must involve a large 
component of strain, which is the secondary and not the 
primary order parameter of the structural transition. 
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